Density functional formalism in the canonical ensemble 
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Density functional theory, when applied to systems with T 7^ 0, is based on the grand canon- 
ical extension of the Hohenberg-Kohn-Sham theorem due to Mermin (HKSM theorem). While a 
straightforward canonical ensemble generalization fails, work in nanopore systems could certainly 
benefit from such extension. We show that, if the asymptotic behaviour of the canonical distribution 
functions is taken into account, the HKSM theorem can be extended to the canonical ensemble. We 
generate A''-modified correlation and distribution functions hierarchies and prove that, if they are 
employed, either a modified external field or the density profiles can be indistinctly used as inde- 
pendent variables. We also write down the A^-modified free energy functional and prove that its 
minimum is reached when the equilibrium values of the new hierarchy are used. This completes the 
extension of the HKSM theorem. 
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Density functional theory (DFT) is, undoubtedly, one 
of the more reliable and established tools in condensed 
matter theory. It has successfully been used in a great 
variety of classical systems II) as, e.g., uniform and non 
uniform systems in simple and general Q| fluids, 

confined fluids |^,^, melting and freezing |^], interfaces 
1^, etc. as well as in the calculation of electronic prop- 
erties in all sort of systems . This impressive amount 
of work directly descends from the pioneering work of 
Hohenberg, Kohn and Sham and its extension to 
non-zero temperature for systems described in the grand 
canonical ensemble (GCE) by Mermin ijllj. Roughly 
speaking, the Hohenberg-Kohn-Sham-Mermin (HKSM) 
theorem states that, either the external potential or the 
density profile can indistinctly be used as independent 
variables and that the thermodynamic grand potential 
reaches its minimum value when the equilibrium density 
profile is used. Therefore, it is in the foundation of all 
sort of variational principles. The failure to implement 
a straightforward canonical ensemble (CE) extension is 
already well known, referred to in Q , and can be traced 
back to the fixed N constraint (see eq. @))- As the 
GCE is an extension of the CE, an obvious question is 
where the need for such an implementation exists. In 
relation to that question, we can mention that experi- 
ments done in porous glasses (with a mean pore radius 
r « 20-301 ) fill, simulations |l|l and DFT studies | 



show the interest of having a mesoscopic DFT in the CE 
in order to study the statistical mechanics of finite closed 
systems like fluids in spherical pores Q. In this letter 
we consider a general mixture of p components (labeled 
by greek indices) with Na particles in a box of volume 
V and under the influence of one body, species depen- 
dent, external potentials V"i^''(x). Following Ramshaw's 
work we will show how, by stripping the canonical 
correlation functions off of its asymptotic behaviour, an 
N-modified set of distribution and correlation function 
hierarchies (which we also show can be defined by func- 
tional derivatives of the canonical partition function Qjq 
with respect to a modified external potential) is intro- 
duced and it paves the way to prove the HKSM theorem 
in the CE. More specifically, we do that by introducing 
an iV-modified free energy which is minimized by the A'^- 
modified density profiles. Here we present the essential 
results with the pair correlation function case as a pro- 
totype and a complete derivation will be published in a 
paper in preparation. 

As usual the distribution functions and full 

distribution functions, the truncated correlation func- 
tions associated to them and direct correlation functions 



hierarchies 
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(c^''({x}) respectively) can be defined by functional 
derivatives as 
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where a — (ai, . . . , cts), {x} = (xi,...,Xs) are the 
species and coordinate sets, /? = {kT)~^, F^^'^ is the 
excess free energy and $q(x) = — /3vi^^(x). The GCE 
definition is formally the same with the grand partition 
function S replacing Qjy. We will consider that the full 
distribution and correlation functions n^^'^i^*) are the 
fundamental quantities (together with c'^*'') in the sense 
that they not only pinpoint the problems in the CE, is 
through these functioncs that we solve them and no new 
information can be gotten from the other functions. It is 
obvious that = fia'^ = ta"^ — £q and that 

<5n(i)(x) [ iiA(x,y)'5'i>A(y)dy (6) 
/ 4A(x,y)c«y-| (1) (1) 



The link with the more usual notation is r^l{x,y) = 
^ali'^^y) + "a'(x)(5„A'5(x-y) where ii^^(x,y) = 
nQ^(x)n^^^(y)/i[j^_^(x, y) and, in this case, the indepen- 
dent variables are the external fields. If the density pro- 
files {ria ^} can also be used as independent variables, 
then an inverse kernel \z,x) = (5$^(z)/(5nL^^(x) 
must exist such that 

^ / £(i'-i(z,x)ti^i(x,y)dx = <5,,5(z-y) (7) 

Using the definition of the pair distribution function in 
both ensembles , it can be obtained that 

CE 



Eqs. (0) and (^) are obviously incompatible in the 
CE while there is no incompatibility in the GCE; be- 
cause of the fluctuations, the nucleus t , is invertible. 
Therefore, the Ornstein-Zernike (OZ) equation (obtained 
when c^a(z,x) is defined by i[fa ^(z,x) — — c^q(z,x) + 
SriaS{z — x)/nQ''(x)) is, mathematically speaking, unde- 
fined in the CE and, within that framework, no rigorous 
DFT is possible. In other words, any variational principle 
formulated in the CE with the canonical functions hier- 
archies and having the density profiles as independent 
variables does not have a mathematically sound founda- 
tion. Specifically, this is the problem alluded to in 

The way out we present in this letter is through strip- 
ping the correlation functions f'*' off of its asymptotic be- 
haviour. This stripping defines the new, iV-modified dis- 
tribution and correlation functions, ri^^^ ^'^^^ respectively. 
We then show that f'"-* is invertible and, from i^^'>~^ we 
define A^-modified direct correlation functions c'-'*^ , write 
down OZ equations that link i'^'^'> and c'^^^ and, by proving 
that they can be also obtained by functional diferentia- 
tion of the CE partition function with respect to modified 
external fields $a, we conclude that rj^"), i'*^ c^**^ consti- 
tute a new hierarchical set of distribution and correlation 
functions in the CE without the mathematical restric- 
tions that the more conventional ones have. Lastly, we 
write down a A^-modified free energy in terms of these A'^- 



modified functions and prove that it reaches its minimum 
value when the equilibrium A^-modified density profiles 
{ria ■*} are used. This completes the proof of the exten- 
sion of the HKSM theorem to the CE. Here we will sketch 
the main steps using the pair correlations case as a pro- 
totype and write down the general results. The complete 
derivation will be reported in a paper in preparation. 
The irreducible two bodies behaviour is exclusively 

(2) 

contained in h^j^ and we will analyze its asymptotic be- 
haviour. The conditional probability of finding an a par- 
ticle in X when a A particle is fixed in y can be written, 
when x and y are very far away, as 

n„A;oo(x|y) = nii)(x) + A0,(y) (9) 

Here A0A(y) is an unknown proportionality factor (ob- 
viously {A,y} dependent); it is not the difference in den- 
sity due to anchoring a A particle. Notice also that the 
derivative is with respect to pa, not p\. The reasons are 
not only the particle exchange symmetry, they are also 
related to the fact that one of our goals is to obtain a full 
hierarchy of A'^-modified functions (otherwise, in eq. ( p^ ) 
a modified external field could not be defined). There- 

(2) 

fore, the asymptotic behaviour of h^^ can be written in 
a fully symmetrical way 
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dlnna' (x) dlnn)^' (y) 



dpa 
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Ta\ is an as yet unknown constant. Next, we define an 

~ (2) 

iV-modified correlation function ( as the correlation 
function stripped off of its asymptotic behaviour. More 
precisely, 



/^i'i(x,y) 



l(x,y) 



■ = y 
(11) 



It must be noticed that, as the stripping is done through 
a separation of variables the irreducible two bodies com- 
ponent is not affected, only the long range behaviour (due 
to the fixed N constraint) is isolated. As a consequence, 
we want to emphasize that hi^\ is the correlation function 
with a truly irreducible two-bodies behaviour, not h\y^- 
Also, the excluded volume effects are not altered and the 
exchange symmetry imposes that Y a\ — -Therefore, 
we can define an TV- modified full truncated correlation 
function 



iS(x,y) 



tS(x,y) 
-Ci(x,y) 



J. 31nn;,^)(x)ain»'/'(y) 



and verify that it satisfies 



X = y 
(12) 

(13) 



This shows that is indeed invertible. Also, Tax can 



be written as a functional of h 



(2) 



y2 



dx7iii)(x) 



<5aA+ / dynW(y)/ii'i(x,y) 



(14) 



i.e., a sort of averaged compressibility. This recovers the 
classical Lebowitz-Percus results |lq] in a form suitable 
for our purposes. Replacing eq. (|12|) in (^) and using the 
chemical potential definition, it can be obtained that 



<5n(i)(x) 



E 

A 



iti(x,y)<5$A(y)dy 



where <i>A is the modified external field 
<5$A(y) = (5$A(y) + /3i5/iA 
If we write that, by definition. 
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(5a,,(?(z - x) 



nL'^(x) 



(17) 



the iV-modified functions t 



r(2)- 



therefore, K 



(2) 



-('2') 

^ and cj,a are linked by an OZ equation with 
c},Q playing the role of an A'^-modified direct correlation 
function. 

Now we define an iV-modified intrinsic free energy 
functional by 



1 7(2) 1, 

'*aA obey eq. 



and, 



(18) 



where U, Pn are the kinetic, potential energy and iV- 
bodies probability distribution function in the CE respec- 
tively with the average being a CE one and minimize the 
functional 



(3n 



E 



(x)$„(x)dx+/3F 



(19) 



It is easy to prove, using the Gibbs-Bogoliubov inequality 
JlSf , that the minimum is reached when the equilibrium 
profiles {fSa^} are used. This completes the proof of the 
extension of the HKSM theorem to the CE. It can also 
be seen that the A^-modified direct correlation function 
hierarchy {c^s)} can be started by 



5(1) 



(x) = /3- 



(5ni'^(x) 



And in particular, we have 



^qA 



(x,y) 



5~c^^\-K) _ 5„A<5(x-y) ,5-l„(x) 



fl'a 



(x) 



5n« 



(y) 



Therefore, we have proved that all of our A^-modified 
functions can be generated by functional differentiation 
with the formal structure of eqs. (||), (^ and (||). 

Summarizing, we have proved that the HKSM theorem 
of the GCE can be extended to the CE and, in the process 
of doing that, we have defined a new set of TV-modified 
distribution and correlation functions that can be built 
using the same rules and formal structure of the conven- 
tional canonical functions with the proviso that the set 
of modified external fields must be used. A paper with a 
more detailed derivation together with some applications 
is in preparation. 
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